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Abstract 

A mathematically motivated fitting function for analyzing supernovae distance-redsliift data 
is suggested. It is pointed out that this function can be useful in analyzing SNAP data. To 
illustrate our method, we analyze current experimental data which includes 78 supernovae and 
20 radio galaxies. From these data we obtain vacuum equation of state in a model independent 
way. We show that the flat Universe model with a cosmological constant is in a good agreement 
with observations. 

1 Introduction 

There is now strong evidence from observations pQ, [2], [3] that the Universe is dominated by a 
nearly homogeneous component — the dark energy — which is causing the cosmic expansion to 
accelerate 0], At present, the distance (magnitude)-redshift relation data is the main source of 
information about dark energy. The new improved data about dark energy will be obtained in the 
nearest future from SNAP experiment. The SNAP satellite will observe roughly 2000 supernovae 
(SNe) a year for three years with very precise magnitude measurements out to a redshift of z = 1.7 
The problem is to establish the energy and pressure as a function of time. In this paper we 
will focus on the new method for analyzing experimental data which can be useful in solving the 
problem. 

The formal solution of this problem is well known. First of all, one needs to obtain the distance- 
redshift function r(z) from experimental data. This function is usually referred to as a dimensionless 
coordinate distance function. This procedure was carried out in pQ, [2] for SNe and in for radio 
galaxies (RGs). Then it is necessary to find the first two derivatives of this function: 
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Now, dark energy (vacuum) equation of state can be expressed through experimental data: 
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The total energy density parameter Q and the total matter density Qm are obtained from indepen- 
dent observations. Recent CMB measurements by WMAP give = 1.02±0.02, Q M = 0.27±0.04 g). 

The main problem in this approach is that it requires numerical differentiation of noisy exper- 
imental data r eX p(z). For any set of experimental data this procedure gives r' exp (z), r" xp (z) with 
many non-realistic fluctuations which have large amplitude in comparison with average values. For 
this reason it is necessary to use smoothing of initial experimental data. The simplest smoothing 
procedure is fitting. However, fitting precision is completely determined by precision of the fitting 
function. Several studies have examined different fitting functions. For example, |§] consider the 
polynomial fit of r(z) and |10j the polynomial fit of e vac (z). Three parametric fitting function for 
apparent magnitude m(z) is suggested in ^3]. The authors of J2| proposed specially constructed 
fitting function reflected some properties of r(z). The model dependent approaches were considered 
in [3|, The model independent approach where r(z) was fitted locally with a polynomial fitting 
is considered in 7 . In our opinion, the problem of constructing the fitting function has a unique 
solution. In what follows we suggest physically and mathematically motivated function allowing one 
to perform the model independent analysis of experimental data and to constrain the dark energy 
equation of state. 



2 Motivated fitting function 

We suggest the motivated fitting function (MFF) which has the following features: 

• MFF can be derived from Einstein equations i. e. MFF is not a phenomenological function; 

• MFF is model independent; 

• a regular mathematical procedure transforms MFF into formally exact solution of Einstein 
equations; 

• for any cosmological solution a(z) with z < 2, MFF approximates a(z), a'(z), a"(z) with any 
given precision. 

With this function, one can increase precision of the model independent approximations of r exp (z) 
and two of its derivatives r' (z), r''(z) systematically. As it is turned out, MFF is completely 
determined by the following condition: the only function we need to approximate is scale factor a(t). 
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It is important that no any other approximation is required. Let us write cosmological solution and 
time-redshift relation 
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The integral (JIJ is approximated as: 
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where x z = x(z) is an exact solution of J3J with x — > when z — > 0. The polynomial coefficients 
a n are considered as fitting parameters. It is obvious that when N — > oo the polynomial coefficients 
can be expressed through derivatives of the cosmological solution a n — ► a^™- 1 /n! at the moment of 
observation. In this limit © - © are formally exact equations. Unfortunately, one cannot reach 
this limit in analyzing current experimental data. For example, the first three coefficients approach 
the first three derivatives starting with N — 8. However, the different cosmological solutions are 
approximated with a good precision staring with N = 4. 

One of the great advantages of MFF approach is that the integral can be expressed through 
elementary functions after some redefinition of the fitting parameter set a n : a,2, 03, 04 • ■ ■ — >, a, b, 
c. . . . (The redefinition procedure will be described later.) The only non-trivial point is that the 
fitting procedure must include exact solution of equation J3J. However, one can use the fact that 
the function x z = x(z) = x(z) has the argument z = z/(l + z) < 1. So, there exists an iteration 
procedure for this function which is always reduced with any given precision. 

In spite of the fact that the result of integration JSJ) becomes cumbersome with N — > 00, there are 
no fundamental difficulties in working with MFF in case of arbitrary N. Besides, explicit calculation 
of © is not necessary. The fitting computer program (which is included into Origin 7 ^3) allows 
us to perform numerical calculation of © during fitting procedure. This approach becomes efficient 
starting with N > 5. 

To perform these calculations we need to restrict polynomial order N to a minimal one. This 
restriction is fixed by the following two conditions. The first condition is that the polynomial should 
reproduce a(x) and its two derivatives a, a with respect to x = Ho(t — to). The second is that the 
fitting curve should be stable against increasing the polynomial order. The second condition can be 
satisfied during fitting procedure. To meet the first condition one needs to study the mathematical 
features of the approximation. We can find the vacuum equation of state from Einstein equations 
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only if we have a good approximation for the first and second derivatives of a(x). The parametric 
representation of the equation of state e vac — e vac (a), p vac = p vac (a) is derived by substituting 
x = x(a) into JBJ|. 

We studied polynomial approximations of cosmological solutions for quintessence models e vac = 
wpvac with —1.6 < w < —0.8, £Ide = 0.73 ±0.04. During numerical experiments it has been shown 
that N = 5 MFF approximates any cosmological solution with relative error \5e vac (x) / e vac (x)\ < 
\Sp V ac(x)/p vac (x)\ < 0.01 for z < 2. If we take N = 4 (with 3 fitting parameters), we have 
\8e vac (x)/e vac (x)\ -0.03, \8p vac (x) / p vac (x)\ -0.1. 

High precision and universality of the polynomial approximation of the scale factor and its deriva- 
tives means that this approximation does not represent any model. So, approximated polynomial 
with the coefficients obtained by experimental data fitting can be considered as a model independent 
cosmological solution. 

Let us write approximated function r(z) for N = 4 explicitly. Polynomial coefficients are rede- 
fined in the following way: 
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Equation 10} can be solved through the following iteration procedure 
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where K is the number of iterations with \x z (^x+i) — x z(K) \ < 10 . 
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3 Experimental data fitting 



As we assume, MFF can be useful in analyzing SNAP data. However, in this work we consider 
current data to illustrate suggested MFF method. Namely, we consider 78 SNe data set (Table 3 
from 0) and 20 RGs data set (joint fit in Table 1 from[7j) for z max — 1.79. Some of the SNe were 
observed more than once by independent observations (the corresponding events will have greater 
weights). In this work we do not average out this events, so we consider the array of 112 points. 

The fitting result is shown in Figure^ We used instrumental weighting and Levenberg-Marquardt 
method |14|. The obtained parameters are 

a = 0.46183;°™ , b = 0.53834+.°,™ , 

(9) 

c = -0.47905183}$ X^/DoF = 1.07499+°™ ■ 
Upper and lower values in correspond to parameters for upper and lower curves which bound 
CL68% region. 

This data set has been fitted by MFF with N — 3 and N = 5. It is turned out that fitting curves 
are almost coincide in the region where z < 1. In the region z > 1 MFF with N = 3 does not work 
satisfactorily (for example the corresponding x 2 is greater than those in 10 •) On the contrary in the 
region z > 1.5 MFF with N = 5 is very sensitive to single experimental points and has non-realistic 
wave-like behaviour. For this reason MFF with N = 4 (3-parametric function) is optimal for the 
analysis of the experimental data. Note, that for the SNAP data, optimal N will be greater (most 
likely N = 5). 

One can see from © that we have a set of fitting curves with almost the same \ 2 - Every fitting 
curve corresponds to a cosmological solution and we need to find the true solution. To do this we 
attempted to average experimental data. We considered two types of averaging. In the first type we 
calculated average values summing data on the different small Az regions. In the second type we 
divided z axis into three or four regions, fitted data independently on each region and then combined 
resulting curves using MFF. These methods gave almost the same result: averaged curves lie inside 
CL68% not far from initial fitting curve. After the long numerical experimentations we concluded 
that initial fitting curve is a better candidate for reflecting the true cosmological solution. 



4 Vacuum equation of state 

From the fact that geometrical set of averaged data points lie not far from initial fitting curve it 
follows that MFF with N > 4 can reconstruct structural features of the true solution at the different 
redshift regions. This allows us to obtain the first derivatives r'(z),r"(z) directly from analytical 
expression for MFF. For this reason rf it (z) can be regarded as the mean values and Ar syst (z) — 
\ r exp( z ) ~ r fit{z)\ as the systematic errors. So, each experimental point can be parameterized by the 
statistical confidence parameter 

, _ Ar stat (z) + Ar syst (z) 

r fit( z ) 
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where Ar stnf (z) is the statistical error (denoted in [7] as <j{yi)). We suggest to use this parameter for 
estimating various physical quantities Tfu{z) calculated using fitting curve. Here we can introduce 
conventional errors AJ r (z) = w(z)\Tfit{z)\. 

Let us turn to vacuum equation of state and consider the function 
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which is derived using data fit (see FigureEJ- When z ~ 2 it is necessary to use both functions in (fTTT|) . 
It should be emphasized that functions plotted in Figure were obtained in a model independent 
way. 

Since p e f / (a) changes slowly it is reasonable to suppose that flat Universe model with a cosmo- 
logical constant could be a possible cosmological model. More detailed information can be obtained 
through fitting e e ff(a), p e ff(a) with the following model 
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which has three fitting parameters £Ide, n, f2. As a result we have 

fl DE = 0.64290 ± 0.07385, n = -0.14445 ± 0.23847, 

Q = 0.98801 ±0.10433, 
w = P vac ^ = -1.04815 ± 0.07949, 

£vac Xp') 

X 2 = 0.04584. 

For non-flat LCDM model (n = 0, Qde = ^a) we have: 

Q A = 0.68973 ± 0.01897, fi = 1.05293 ± 0.03068, 
X 2 = 0.04753. 



(12) 
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One can see from (I12|) . (|13|l that these results are in a good agreement with WMAP experimental 
data [Hj • We hope that presented MFF method can provide an independent check to the traditional 
cosmological tests. 
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Figure 1: Dimensionless coordinate distance r(z) as a function of z. The solid line is the best fit. 
The dotted curves indicate CL68%. Experimental data are taken from [7] (full data set). 




Figure 2: Effective energy density and pressure as a function of a obtained through model indepen- 
dent analysis of experimental data using MFF. Fitting curves correspond to the model with 
parameters from l|12l) . 
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